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ABSTRACT 


An expression for obtaining the lifetimes of ar- 
tificial satellites in circular orbits is developed in 
this paper. A complete derivation of the method is 
presented to allow the user to evaluate its assump- 
tions according to specific needs. The accuracy of 
the developed method is verified using Earth and 
Mars as examples and comparing the results to the 
results obtained from numerically integrated and ap- 
proximate analytical trajectories. Characteristic 
altitude histories are presented for the trajectory 
methods and the described analytical solution. Ex- 
pressions for computing lifetimes of elliptical orbits 
are also included, with the necessary graphical pres- 
entations to provide rapid solution of these expres- 
sions. 
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A GENERAL ANALYTICAL METHOD FOR ARTIFICIAL-SATELLITE 


LIFETIME DETERMINATION 

By Frank Garcia, Jr. 
Manned Spacecraft Center 


SUMMARY 


Existing analytical techniques for the determination of satellite lifetimes incor- 
porate assumptions which result in simplified expressions at the expense of generality. 
This paper develops an expression for obtaining lifetimes of circular orbits which car- 
ries terms in the order of the cube of the quotient of the altitude over the radius of the 
planet which is necessary when considering planets smaller than Earth. A complete 
derivation of the method is presented to allow the user to evaluate its assumptions ac- 
cording to his specific needs. The accuracy of the developed method is verified using 
Earth and Mars as examples and comparing, the results to the results obtained from 
numerically integrated and approximate analytical trajectories. Characteristic altitude 
histories are presented for the trajectory methods and the described analytical solu- 
tion. Expressions for computing lifetimes of elliptical orbits are also included, with 
the necessary graphical presentations to provide rapid solution of these expressions. 


INTRODUCTION 


Any artificial satellite which enters the atmosphere of a planet during a portion 
of its orbit loses a small amount of energy per revolution as a result of atmospheric 
drag. This loss of energy may be categorized as occurring during two phases: (1) that 
which takes place when the orbit has appreciable eccentricity and (2) that which takes 
place when the orbit is nearly circular. If the orbit of a satellite has an appreciable 
eccentricity, most of the air drag is encountered in the vicinity of periapsis. Thus, the 
satellite suffers a slight velocity retardation as it passes through periapsis. This re- 
tardation leads to a distance decrease at apoapsis, while scarcely altering the periap- 
sis distance. A good physical example of this phenomenon is observed in Explorer I, 
which has suffered a velocity decrease at periapsis of 14 n. mi. and a distance 
decrease at apoapsis of 695 n. mi. in a period of 9 years. Thus, air drag tends to 
circularize an orbit of initial appreciable eccentricity. An orbit of very low or zero 
initial eccentricity will manifest its energy loss as slight contractions of the orbit. 
These energy losses may be used to assess satellite lifetimes using analytical tech- 
niques enhanced by simplified but realistic assumptions. 



A critical review of currently available mathematical tools to estimate artificial- 
satellite lifetimes is presented in reference 1. The technique presented in this paper 
has been developed from the same energy considerations and differs principally in the 
treatment of the mathematics and the inclusion of higher-order terms. It is considered 
that the resulting expressions now provide a somewhat more accurate expression than 
similar solutions for Earth alone. In particular, the inclusion of higher-order terms 
makes the satellite lifetime expression amenable for use with smaller planets such as 
Mars. 
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SYMBOLS 

o 

satellite reference area, ft 
semimajor axis, ft 
ballistic number, mJC^A 

drag coefficient 

aerodynamic drag, lb 

2 2 

total energy per unit mass, ft /sec -slug 

eccentricity of orbit 

exponential function 

altitude above planet, n. mi. 

modified Bessel function with argument /3x q 

4 

time conversion factor, 3153. 6 x 10 sec/yr 

2 2 

kinetic energy per unit mass, ft /sec -slug 

satellite mass, slugs 

revolutions 

variable in integration of N 
orbital period, sec/ revolution 
mean orbital period, sec/revolution 
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PE 

2 2 

potential energy per unit mass, ft /sec -slug 

R 

radius of the planet, n. mi. 

r 

radial distance from center of planet to satellite 

T 

normalized satellite lifetime, t/ BN 

t 

satellite lifetime, yr or days 

V 

orbital speed, ft/ sec 

V C 

C () 

circular speed at ( ), ft/ sec 

X 

average altitude, ^H & - H^/52 

/3 

density scale height factor, 1/ n. mi. 

M- 

3 2 

universal gravitational constant, ft / sec 

P 

g 

density of air, slugs/ ft 

Subscripts: 

a 

apoapsis 

H 

altitude above planet, n. mi. 

o 

initial or surface value 

P 

periapsis 

() 

quantity evaluated at point in parentheses 
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Earth 


METHOD OF ANALYSIS 


The analysis is made for orbits of initial zero eccentricity and for those with ap- 
preciable eccentricity. The lifetime of a circular orbit is defined as the time that it 
takes an artificial satellite to decay from its initial altitude to an altitude correspond- 

_ JQ g 

ing to a density of about 1 x 10” slugs/ ft . However, the lifetime of a satellite in an 

orbit of appreciable eccentricity includes the time that it takes it to decay to a circular 
orbit from its initial eccentricity. 
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It is assumed that the gravitational forces present are those caused by a purely 
inverse- square potential. The nongravitational forces are assumed to be caused by a 
nonrotating exponential atmosphere. 


Pjj = P 0 exp(-/3H) 


D ' 2 P H C D AV C 


( 1 ) 


These assumptions are standard in the approach to an analytical expression for satellite 
lifetimes. 

Presented in reference 1 is an expression similar to the one presented in this 
paper for the lifetime in revolutions. Also presented in reference 1 is a similar ex- 
pression which can be obtained from equation (A14) in appendix A by neglecting terms 

2 

in the order of (H/R) and assuming that H » 1//3. A lifetime expression derived 
by using an arithmetic mean is given in reference 2. 
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The method presented in this paper carries terms in the order of (H/R) and 
uses an integrated mean which would reduce error introduced by taking large altitude 
intervals when using the criteria already suggested in the form of an exponential atmos- 
phere. The choice of an integrated mean instead of an arithmetic mean for the 
period P was suggested by the nonlinearity of P with respect to r. Terms of higher 
order were carried to make the lifetime expression amenable for use with smaller 
planets such as Mars. 


Circular and Near- Circular Orbits 

At each point on its orbit, a satellite is subject to drag in the opposite direction 
of the velocity V^,. The circular speed is 



During one revolution, the orbit remains circular with a radius of r if it is assumed 
that the variation in Ar is small. This assumption can be verified by integrated tra- 
jectories which show, if the orbit is initially circular or near circular, that during life- 
time of the orbit this near- circular condition is maintained with eccentricities on the 
order of 0. 0001 or less (fig. 1). Then, for one revolution the change in total energy is 
caused by drag (fig. 2), and this change is commonly expressed as 

AE = / D r d6 (3) 
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Using this energy approach, in appendix A it is shown that the satellite lifetime may be 
expressed as 



2 

where T = t/BN yr-ft /slug. In computing lifetimes using equation (4), it should be 
noticed that 0 and P q are assumed to remain constant between Hj and Hg. A semi- 
log plot of density versus altitude will dictate which and Hg to choose for the best 
results; that is, and Hg should be chosen so that the density scale height is es- 
sentially constant within this interval. The method may be continued analytically in 
this way until the final desired altitude is reached. However, for the sake of expedi- 
ency, equation (4) was set up on a digital computer, with a subroutine which computes 
P q and 13 as a function of altitude. 


Orbits of Appreciable Eccentricity 

A satellite on an orbit of appreciable eccentricity (/3 x q > 2) suffers its greatest 

loss of energy at periapsis, and this loss will have more effect on the decrease of ve- 
locity at periapsis than on the radius there. The net effect would be similar to an im- 
pulse applied at periapsis. Thus, it is assumed that the radial distance at periapsis 
remains nearly constant (fig. 3), making use of the exponential atmosphere 

P H =P 0 ex Pp( r p o - r 
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The restriction on /3x q is caused by the expansion of the modified Bessel function 
I l(/3 x o ), which appears in the expression for lifetime. In reference 3 it is shown that 
the lifetime may be expressed as 


T = 



Kpft) 

y o 



( 6 ) 


where (3 is computed at the initial perigee point p Q and remains constant. Here 

2 

again, T = t/BN yr-ft /slug, and it is shown in appendix B that 



2 

where only terms in the order of e Q are retained. The assumption that the radial 

distance at periapsis remains nearly constant is valid in the physical sense as long as 
BN does not approach zero. Conversely, as BN approaches large values, this as- 
sumption agrees well with the actual situation. 

The parameters f ^/3 x q , e Q ^ and g^(3x^ (shown in figs. 4 and 5) may be used to 
compute lifetimes of satellites in elliptical orbits (ref. 3). 


Sample Computation of the Lifetime of an Earth Satellite 

The lifetime of a satellite in an initial elliptical orbit may be divided into two 
phases. The first phase is that time which the satellite theoretically spends in produc- 
ing a nearly circular orbit (fig. 3), and the second phase is the time that the satellite 
spends in the decay spiral (fig. 1). The sum of the time spent in these two phases be- 
comes the total lifetime. The following is an example to clarify the computation 
method. 


J 
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The planetary constants used are 


|jl = 1. 4020066 X 10 16 ft 3 /sec 

R = 0. 209020 X 10 8 ft 

V n = 25 898. 33 ft/sec 
'-'o® 

B = 10 slugs /ft 3 

H = 2000 n. mi. 
a 

H = 300 n. mi. 

P 

/3 = 0. 0221 8/n. mi. 

p = 2 x 10" 15 slug/ft 3 
p o 



then the necessary constants are 


e =0.185 \ 

o 

/3x q = 18.85 


a = 4588. 7 n. mi. 
o 


= 20591.3 ft/sec 

% 


K | ip 8 = 3225. 7 lb-sec/yr-ft" 

( P o) 


V r Map Z 3 = 6. 3835 yr-ft / slug 

\ ° 


Entering e Q and /3x Q into figure 4 


f (^o’ e o)" 110 


therefore, the lifetime in this phase is 


T x = 702.185 yr-ft / slug 


With H = 300 n. mi. in figure 6 


Tg — 2700 days-ft /slug = 7.397 yr-ft /slug 
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The total lifetime T is then the sum of and Tg, or 

T 1 + T 2 = T = 709. 582 yr-ft 2 / slug (14) 

RESULTS AND DISCUSSION 


Comparison of Results 

The results obtained from the analytical expression developed in appendix A for 
computing circular orbit lifetimes are presented for both Earth and Mars. The results 
for Earth are compared to two trajectory simulations (fig. 6), whereas the Mars re- 
sults are compared to only one simulation (fig. 7). 

The two trajectory simulations are termed as the analytical ephemeris generator 
(AEG) and the numerically integrated point mass program (GEM). The AEG program 
makes use of Lagrange’s planetary equations with an oblate Earth and drag effects. 

The 1962 U. S. standard atmosphere (ref. 4) shown on figure 8 is used, and the oblate- 
ness is described by the second, third, and fourth harmonics of the Earth. The main 
use of this program is not for the prediction of orbital decay; but because it is a rela- 
tively fast program, it may be used for this purpose. The GEM program (ref. 5) con- 
siders the point mass in three-dimensional space and integrates with a fourth- order 
Runge-Kutta integration routine. This simulation also uses the 1962 U. S. standard 
atmosphere. 


Analysis of Results 

As stated previously, the definition of a lifetime as used in this report is the time 
from the initial orbital- insertion altitude to the arrival at that altitude where the den- 

-10 3 

sity reaches about 10 slugs/ft . This value corresponds to about 50 n. mi. for 
Earth and 90 n. mi. for Mars with the models used. It was also noted in the GEM that 
at about this point the magnitude of the flight-path angle increased greatly, and thus, 
the time to impact became negligible as compared to the overall lifetime of the satel- 
lite. 


The lifetime of a satellite in orbit about the Earth, as calculated by the three 
prediction methods using the U. S. standard atmosphere (fig. 8), is presented in fig- 
ure 6 for comparison purposes. The maximum deviation of the analytical technique 
from the numerically integrated trajectory is less than 13 percent. The maximum de- 
viation from the Lagrangian planetary trajectory is less than 25 percent. These devi- 
ations occur at the higher altitudes and decrease rapidly thereafter. At 200 n. mi, 
the deviation of the presented technique from the GEM and AEG is less than 9 and 
19 percent, respectively. 
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The lifetime for a satellite of Mars is shown in figure 7. The results were ob- 
tained using the GEM and the analytical technique with the atmospheric model of fig- 
ure 9, which was obtained from reference 7. Here, the maximum deviation is less 
than 10 percent, proving that the analytical technique can be used to predict the lifetime 
with reasonable accuracy. 


CONCLUDING REMARKS 


Any simulation which requires numerical integration is not only costly but subject 
to accumulative errors when used continuously for a long time. Of the three methods 
used in this report for orbital-decay determination, the closed-form method is by far 
the least expensive, with the numerically integrated point mass program being prohibi- 
tive whenever large lifetimes are expected. The analytical ephemeris generator is ad- 
equate for lifetimes of 2 years or less. 

The deviations of the analytical technique from either the numerically integrated 
point mass program or the analytical ephemeris generator trajectory are within the 
present capability to predict either the atmosphere or the drag coefficient of a body in 
rarefied flow. Therefore, the relatively small deviation from the numerically inte- 
grated point mass program implies that this analytical technique is adequate for orbital 
lifetime prediction. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, October 18, 1967 
981-89-00-00-72 
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Circular insertion 



Figure 1. - Theoretical decay of an initially circular orbit 
perturbed by a planet atmosphere. 
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Figure 6. - Lifetime of a satellite in an initially 
circular orbit about Earth. 


Satellite lifetime, T, years/BN 




APPENDIX A 


CIRCULAR-ORBIT LIFETIME 
For one revolution, the change in energy is 

AE =f Dr d0 = -f | pC D AV c 2 r d 9 (Al) 


or 



h c d a / 


p dd 


(A2) 


12 2 

where D = - ^ pC D AV^ and = p/r. The closed-path integral is taken along a 

circle of radius r with 6 as the argument being swept out (fig. 2). Where an expo- 
nential atmosphere is assumed, these equations become 


1 

AE = - g p.p o C D A exp(-/3H) dd, 


(A3) 


and assuming that during one revolution the change in altitude is small, the loss of 
energy per revolution is 


AE = -pp o C D A7r exp(-/3H) 
The total energy per unit mass is 


(A4) 


E = 


KE + PE = 



J± 

r 


(A5) 


E = - 


2r 


(A6) 
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K the mass of the satellite is m, then 


E = - 


mp 

2r 


(A7) 


and 


AE = mi (i?) = -fis ir 


(A8) 


assuming that m remains constant. 

Since the change in energy is equal to the work done, substituting equation (A4) 
for AE in equation (A8) yields 


Ar = 


2AE r 2 _ 2 tjt 2 
m pl = " 7BNj P o 


exp(-/3H) 


(A9) 


Because - Ar is small, the first variation may be taken as the derivative, and recalling 
that Ar represents the change in radius for one revolution, equation (A9) may be re- 
written as 


9 0-7T 

(R + H) exp(/3H) dH = - p Q dn 


(A10) 


Expanding 


(R + H) 2 


in powers of (H/R) 


(R + H)" 2 




(All) 
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Retaining terms in the order of (H/R; only and substituting equation (All) into 
equation (A10) 


MIMS)*- 4 ®' 

\ 

■ 


exp(fiH) dH 

j 


therefore 


r ‘ 2 [ 1 - 2 ( h ) + ii) - <|) 3 ]^p(OT) dH = d„ 


(A13) 


Integrating the left-hand side of equation (A13) from Hg to and the right-hand side 
from zero to N 


ijexp^Hj) 1 - ( H ! - j) (h + ^ + ^3) + (tt) ( 3 + - 4 (tt) 


-^H 2 ) 


2 6 24 \ ( 2 

r+ r + u 


( 3 + m) ‘ 


R / 75MJ p o N 

(A14) 


where N is the number of revolutions to decay from to Hg. It is now necessary 
to obtain an average period to get an expression for lifetime in terms of time. 


NP = t 


(A15) 


Making use of the first mean- value theorem for integrals, 


2 n / * „3/2 ^ 4 tt / 1 


2 2 

r l_ r 2 


(A16) 
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The integrated mean of the period from ^ to r 2 is P as expressed in a general 

manner in reference 6. The lifetime expression may be obtained by solving equa- 
tion (A15) for N, substituting for P from equation (A16), and substituting for N in 
equation (A14) to obtain the final desired expression 
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APPENDIX B 


ELLIPTICAL-ORBIT LIFETIME 


Reference 3 shows that the satellite lifetime is 


e 2 ( lie 29e 2 _ 

. _ o [ 1 o o 7 

2CK \ 6 + 16 + 8/3a 


Vtt 

c * ‘ e o) 


_ (BN ><Vo exp (' 3x o) exp e Q lle 0 29e Q , 

■ 4 ^ ( p o) W x o) K V 6 + 16 + 8 »< 


which may be simplified by using the following relationships. 


P = 2? rJ— 

o ^ p 


x = a e 
o o o 


V* 


fro C ( a o) 


42e exp e 


X i^xNk • p> 7 S A 1 48 exp e o - 88e o exp e o + 87e o exp e o + w 

' ' ( p o) \ c 
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or 


T = 


KpP 






(^v 



(B8) 


where 


g 



Dx 0 exp(,.lx o ) 

2 W 


(B9) 


f(/3x o , 



K^o) 

48 


48 exp e Q - 88e Q exp e Q + 87e Q 


exp e 

o 


42e exp e A 
o o \ 

+ 


(BIO) 
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